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■ Abstract. In addition to lines originating in a soft phase at ~ 0.8 keV and to cold molecular clouds, the X-ray 

spectra from the Galactic center region also exhibit properties similar to those of a diffuse, thin, very hot plasma 
at 8 keV on a scale of hundreds of parsecs. This phase is surprising for more than one reason. First, such a hot 
QQ ' plasma should not be bound to the Galactic plane and the power needed to sustain the escaping matter would be 

higher then any known source. Second, there is no known mechanism able to heat the plasma to more than a few 
keV. Recently we have suggested that, hydrogen having escaped, the hot plasma could be a helium plasma, heavy 
enough to be gravitationally confined. In this case, the required power is much more reasonable. We present here 
^ I a possible heating mechanism which taps the gravitational energy of the molecular clouds. We note that the 8 

keV plasma is highly viscous and we show how viscous friction of molecular clouds flowing within the hot phase 
can dissipate energy in the gas and heat it. We detail the MHD wake of a spherical cloud by considering the 
different MHD waves the cloud can excite. We find that most of the energy is dissipated by the damping of 
^— ' Alfvenic perturbations in two possible manners, namely by non-linear effects and by a large scale curvature of 

' the field lines. We find that the total dissipation rate depends on the field strength. For fields B < 200/^G both 

, mechanisms produce power comparable to or higher than the radiative losses; for strong fields B > 1 mG, only the 

curvature damping can balance the X-ray emission and requires a radius of curvature Rc < 100 pc; whereas for 
intermediate fields, the total dissipation is more than one order of magnitude smaller, requiring a higher accretion 
rate. We note that the plasma parameters may be optimal to make the dissipation most efficient, suggesting a 
O ' self-regulation mechanism. The loss of kinetic and gravitational energy also causes accretion of the clouds and 

^ , may have significant action on the gas dynamics in this region between the large scale, bar dominated fiow and 

C/3 . the central accretion to the massive black hole. 

a 
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200 it iMuno et alJ l2004j) (see however iRevnivtsev et al 



astro-p h/0510fl5fl for a different vue). 



1. Introduction tions have shown that no more than 10 % of the emis- 
sion c an be attributed to point sourc es l)Ebisawa et al 

The optical emission from the Galactic center is strongly ' ' ' ' ' ' 

obscured by dust absorption. Radio, IR and X-ray astron- 
omy have thus become useful tools to probe this region. 
The present picture of the center region is very complex, xhe Galactic ridge X-ray emission may thus be truly 
showing supernova remnants, H2 regions, star clusters, ex- diffuse and its spectrum is a powerful diagnostic to under- 
tended filaments, magnetism, star forming regions... For gtand its nature. In addition to the bremstrahlung emis- 
more than 20 years, many high energy observations have gion, an intense highly ionized He-like F e emission line 
reported an intense X-ray emission from the Galactic at 6 7 keV was first reported by Ginga llKovama et alJ 
center region (|Worrall et alj [198^ [Warwick et ^ ll985|) . . The continuum and line are observed to extend out 



Although a cosmic contribution has been identified as re- to several kpc (jYamauchi et al. 199ai with a typical scale 

suiting from discrete sources, the remaining Galactic ridge height of less than 100 pc, b ut a strong enhanceme nt is ob- 

emission has not been resolved yet. Some spectral prop- served in the inner 300 pc jYamauchi et alJll990|) . These 

erties are very similar to those of Galactic X-ray point observations have then been confirmed by many others 

sources (|Wang et al.||200a) but recent Chandra observa- and new H- and Hc-likc lines from very ionized elements 

have been resolved with ASCA in the .5 — 10 keV band 

* e-mail: belmontOcea. fr ijkaneda et al...l997.) . The study of these Mg, Si, S, Ar, 
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Ca and Fe lines indicates that the spectrum cannot be 
modeled with a single-temperature thermal plasma. As a 
result two classes of models have been proposed. 

On the one hand, it has been proposed that 
the spectrum could result from a single-temperature 
plasma assoc iated wi th non- or quasi-thermal mechanisms 
ilValinia et al.i.200a:.Tanaka et al.ll2000HMasai et all2002t 



00:. 

iDoeiel et al.l2002allJ ). However, besides the intrinsic prob 
lems of the different models, recent observations with the 
Chandra satellite have pr ovided spectra mo re consistent 
with a thermal origin (see iMuno et al.ll200^ for a review 
of these models and a comparison with results of Chandra 
observations). 

On the other hand, authors have interpreted the spec- 
trum as originating from two different spe ctral compo- 
nents jKaneda et allllQQTtrMuno et al.ll2004 : a soft com- 
ponent at kT ~ 0.8 keV and a hot component at kT ~ 
8 keV. The soft plasma properties are compatible with a 
supernova origin: its spatial distribution is patchy and su- 
pernova shock- waves are known to generate temperatures 
of ~ 1 keV. 

The origin of the more diffuse 8 keV plasma is puzzling 
and raises several questions. Assuming that the plasma 
has solar abundances, i.e. is essentially made of hydro- 
gen, it has been noted that its temperature is too high 
for it to be grayitationallY bound to th e Galactic plane 
ijKovama et al.lll996l: iMuno et alJ 12004 - The power re- 
quired to compete with the energy losses associated with 
this escaping matter exceeds the power of any known 
source. An average supernova rate 10^ or 10^ higher than 
in the rest of the Galaxy would for instance be required 
to heat the ho t plasma before it escapes the central region 
(derived from lMuno et al.ll2004|) . 

More recently, we have noticed that, at the inferred 
temperature and density of the hot phase, hydrogen 
ions are weakly collisional with other ions so that they 
can leave the Galactic plane without dragging other ele- 
ments with them, leaving a he lium plasma that is heav- 
ier and confined by gravity ijBelmont et al.l l2005|) . As 
this plasma does not escape, energy losses are dominated 
by radiatio n which occurs on much longer time scales: 
~ 10^ yr (iM uno et al."20i}J). From the Ginga data of 
lYamauchi et al,. L1990i|) we derive a peak emissivity of 
1.4-4.6x10'^^ erg s~^arcmin^^ in the central region. This 
result is rather consistent with recent observations by 
Chandra of the inner 20 pc which have rep orted a local 



lumin osity of L=5.-9.x 10'^^ erg s ^arcmin^ I M^mo et alJ 



llo^). The Ginga value is an average over a large area 
whereas the Chandra field is more central and much 
smaller. This could explain the slight difference. With the 
intermediate value, the total emission in the whole central 
region is L « 4. x 10'^^ erg s^^. This power can be pro- 
vided by reasonable sources. Nevertheless, a mechanism 
must still be found that can heat the gas up to 8 keV. 
Supernovae for example could provide enough energy to 
balance the radiative cooling, but, although their tempera- 
ture depends on their age and the external pressure, SNRs 
have not been observed to temperature higher than ~ 1-3 



keV after a few hundred years, significantly smaller than 
the required 8 keV. 

In this paper, we assume that the 8 keV emission in- 
deed originates in a diffuse helium plasma in the Galactic 
center and we investigate the idea that the heating in this 
region could be provided by the friction of cold clouds with 
the surrounding gas. This idea is based on two main facts. 

First, many molecular clouds have been observed in 
CO, CS, NH3 in the central region. The dense H2 clouds 
of about 10 pc size are observed to form in a ring at about 
180 pc from Sgr A*, then detach and spiral inward with 
a significant velocity relative to the surrounding medium, 
typically 100 km s~^. The large number of moving clouds 
represents a huge reservoir of gravitational and kinetic en- 
ergy that may be used to heat the plasma. 

Second, plasmas at temperatures as high as 8 keV can 
be highly viscous ( Braginskii, 1 965,) . This viscosity is how- 
ever different from the usual one characteristic of neutral 
gases. Indeed, in the conditions of the Galactic center, any 
reasonable magnetic field fully inhibits the usual shear vis- 
cosity and so reduces by orders of magnitude the efficiency 
of the corresponding viscous dissipation. However, we will 
show that the remaining bulk viscosity, which acts on com- 
pressional motions, is sufficient to apply a significant net 
viscous stress on the moving clouds. The associated power 
is dissipated in the plasma and heats it. 

The present paper is organized as follows: first, we de- 
scribe in sect. |21 the Galactic center arena, its molecular 
and magnetic content; then, we present the main charac- 
teristics of the viscosity acting in the Galactic center in 
sect. 131 we describe how viscous friction with cold molecu- 
lar clouds can dissipate part of the cloud's kinetic energy 
in sect. 01 and last, we discuss the efhciency of this heating 
in providing the energy required to heat the medium in 
sect. O 



2. The Central Molecular Zone 



The Central Molecular Zone (CMZ, iMorris fc Serabvnl 
covers the first hundreds of parsecs and differs from 
the outer part of the Galaxy by many features. Diffuse 
radio, IR and X-ray e mission are for exam p le higher there 
than anywhere else liBroean et all l2003t iBennett et alJ 
*1994t lYamauchi et al.ll99Cll) . The molecular content of this 
region is higher than in the outer parts of the Galaxy 
and about 10% of the Galactic molecular gas in con- 
densed in the CMZ. The outer limit is composed by a 
dense torus of gas at about 180-200 pc. At this radius, 
the gas density drops from 200 M0pc~^ to 5 Mqpc~^. 
Last but not least, many non-thermal filaments are ob- 
served in the first 300 pc, which are not observed any- 
where else in the Galaxy. In this paper, we focus on a 
region of 150 pc radius and 70 pc high below and above 
the disk, which r epresents the main volum e of the Central 
Molecular Zone l|Morris Sz Serabvnll99i but does not in- 
clude the dense torus. We will refer to the 'inner CMZ' for 
this region, which corresponds roughly to the region where 
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the filaments are observed. We assume that the distance 
to Sgr A* is 8.5 kpc. 

Our model relies on tapping some of the kinetic en- 
ergy of the numerous molecular clouds in the CMZ. 
The statistical properties (size, mass, internal velocity 
dispersion...) of these H2 cloud s have been studied in 
many details llBally et all llQSSt lOka et all Hool l200lt 
iMivazaki fc Tsuboil hOO^ . These surveys reveal a very 
complex medium whose nature, structure and kinematics 
is not fully understood yet. For instance, their formation 
and binding processes are still uncertain. 

More than 150 clouds have been identified but, because 
of projection effects, several clumps can be observed as one 
single clou d, so that the actual cloud number cou l d be a 
bit higher l|Ballv et alJll988l) . iMivazaki fc Tsiiboil l|200nl) 
reported 15 9 clouds a t Gala ctic longitudes — 1° < Z < 
+ 1''42' and lOka et all l|200l|) detected 165 clouds in the 
region — 0°.8 < I < 1°.7. However, several of them could 
belong to the 180-pc molecular torus a nd it is not clear 
howmany really belong to the inner part. ILaunhardt et alJ 
estimated a volume filling factor of a few percents 
for this cold phase in the region r < 120 pc, \z\ < 50 pc. 
Assuming the clouds have a 10 pc size, we derive a corre- 
sponding number of 13 5 clouds This region considered by 
ILaunhardt et alJ l(2002l) is smaller than the region we are 
interested in but again, several of these clouds could be- 
long to the front part of the 180-pc molecular ring on the 
line of sight. Finally, we consider that at least 100 clouds 
belong to the inner CMZ. 

The cloud size is fairly well constrained even if there 
is uncertainty in the cloud size definition. Clouds are ob- 
served to have radii between 1 pc and about 10 pc. Some 
clouds of a few tens of parsec size are also observed but 
they could be comp lexes of se veral smaller den se clouds. 
IMivazaki fc Tsiiboil (jioOO) and lOka et all l|200ll) derived a 
mean radius of 3.7 pc and 6 pc respectively, which leads to 
use 5 pc as a typical radius. However, clouds are thought 
to have very complex shapes and structures. They are of- 
ten very extended with arc-, shell- or lane- shapes. Also, 
many large clouds are composed of several smaller clumps. 
Gas is observed to fill the space between the clumps with 
densities < 10'^ cm~'^, lower than the mean density of 
the clumps (10'' cm~^). This gas is likely to give them 
a common dynamics. In that sense, the clouds may have 
a fractal structure sim ilar to that of the standard ISM 
i Falerarone et alJll99lh . 

The cloud velocity is one of the less constrained param- 
eters. Careful studies of the velocity distribution suggests 
that most of the molecula r gas is in a torus at about 180 pc 
from the Galactic c e nter llKaifu et al.lll972l:IScovillelll972l: 
iBinnev et J] Il99ll iMorris fc Serabvnl Il99(t . Molecular 
clouds of about 10 pc size are thought to form in this 
torus, then detach and spiral inward with a significant ve- 
locity relative to the surrounding medium. The observed 
velocities are in a r ange between Vr, sR— -250 km/s and 
Vlsr= +250 km/s ijBallv et alJll988() . The clouds have a 
huge cloud-cloud velocity dispersion so that it is difficult 
to get a reliable rotation curve. The cloud velocity dis- 



persion within giant molecular complexes is about 30-50 
km/s. This figure has been directly measured in various 
surveys and is in good agreement wi th the observed scale 
height of 50 pc of these complexes ijOka et al.lll998l) . As 
the hot gas is expected to have very smooth properties, 
this velocity dispersion is a minimal value for the cloud ve- 
locity relative to that gas, whatever the rotation profile is. 
Besides this dispersion, the rotation profile is not consis- 
tent with uniform circular rotation. iBinnev et al.l l(l99lh 
have worked out that the gas dynamics of this Central re- 
gion is governed by a bar-potential. They suggested that 
the molecular ring at 180 pc is a transition region, cor- 
responding to the Inner Lindblad Resonance of the bar: 
beyond this torus, the gas orbits along ellipses aligned 
with the bar axis (the so called Xi orbits) and inside this 
torus, molecular clouds orbit on very elongated ellipses 
perpendicular to the bar axis (the X2 orbits). Because of 
its high temperature, the 8 keV plasma is not likely to 
respond much to the bar potential. Moreover, if the per- 
vasive magnetic field is generated by azimuthal currents, 
they must be localized in the molecular ring, providing 
the gas pressure gradient necessary to balance the mag- 
netic pressure gradient. One can thus expect the magnetic 
field lines to rotate at the torus speed and drag the hot 
plasma with them. As a result, the relative velocity be- 
tween the clouds and the surrounding material can be ex- 
pected to be a substantial fraction of the cloud orbital 
speed (typically 200 km/s). Furthermore many clouds are 
observed with a forbidden velocity, i.e. velocity of sign 
opposite to the rotation direction {I < 0, Vi^sR > and 
I > 0, y[,sR < 0) and t hese velocities can be as high as 130 
km/s ()Oka et al.ll998l) . which would imply a huge velocity 
relative to the ambient field. And last, the interpretation 
of filaments as resulting from the interaction of moving 
clouds with a pervasive magnetic field (as discussed be- 
low) requi res a relatiye velocity of 50-150 km/s, typically 
100 km/s (iBallv et al .111 9881: Is^ie 20051 Will will use 100 
km s~^ as the mean cloud velocity. 

Many vertical filamen ts have been reported in 
the Galactic center region llYusef-Zadeh fc Morris! Il987l: 
iLaRflsa et al.l l200nl) . These filaments are non-thermal 
and magnetized, which indicates that the central re- 
gion is magnetized too. The strength and topology of 
the magnetic field is a debated issue. On the one hand, 
it has been suggested that the filaments may trace 
a strong (B ~ '^G) ve rtical pervasive magnetic field 
((Morris fc Serabvnl Il99fil) . Because in most cases they 
are associated with molecular clouds, these extended 
structures have been suggested to result from t he inter- 



'l98^; 


'Morris & Yusef-Zadeh" 


'l989; 


'Lesch fc ReichI Il992l: 


.Rosso & Pelletier 1993; Serabyn & Morris 19941. To drive 



vertical currents or guide accelerated particles along the 
filaments, the field has to be vertical. Its strength has been 
deduced from the observation that the field lines are not 
distorted by the cloud motion. This observation implies 
that this motion is sub-Alfvcnic. Previous estimates as- 
sumed that the magnetic perturbations propagate in a 
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medium of density comparable to that of the clouds (n 
~ 10^'' cm^'^) and estimated mG strengths. However, if 
the Alfven waves propagate in the faint hot diffuse plasma 
(n ~ .1 cm~3), then, to keep the filaments straight over 
a height of 5 times the cloud's radius (wa > St^c), this 
only requires: B > 100/iG. A strong magnetic field is also 
invoked to confine the clouds. Indeed, the latter are not 
massive enough to be confined by gravity. They need an 
external pressure to balance the internal turbulent veloc- 
ity. The hot plasma pressure is not strong enough to bal- 
ance the ram pressure of the clouds but a pervasive field 
with B ~ 0.5 mG could do the job ijMivazaki fc Tsuboil 
|2003). On the other hand, it has also been suggested that 
the non thermal filaments may not trace a strong field but 
may rather correspond to local enhancements of a mean 
field in pressure e quilibrium. The field could thus be lower 
than B ^ 1 mG (|LaRosa et ahlEoOOl) . Recent considera- 
tions on the diffuse non-thermal radio emission gave upper 
limits for the mean field. If these results are confirmed the 
field strength must be lower than ^ lOO/iG I LaRosa et all 
|200,'t(i . In a hot helium plasma at 8 keV, the equipartition 
field is about B ~ 100//G. In this paper, we use this value 
as a reference but we also discuss the consequences of 10 
/iG-1 mG fields. 



3. The Braginskii viscosity 

We review in this section the viscosity of magnetized plas- 
mas, which is very different from the hydrodynamical one, 
usual in neutral gases, and needs to be discussed. 



3.1. General properties 

Let us first consider a non magnetized plasma. Its viscosity 
is isotropic and it is thus governed by the sole coefficient : 



rjo = -96 



5/2 



(1) 



where rrii and eZi are the mass and ch arge of ions consti- 
tutive of the plasma l|Braginskiilll965|) . We note that the 
temperature dependance is stronger than for neutral gases 
(oc because the Coulomb cross section depends on 

the particle velocity whereas hard sphere- type collisions do 
not. This implies that a 8 keV helium -'^ plasma is highly 



viscous: 



77 = 630 g cm s 



8 kcV 



5/2 



(2) 



The corresponding kinematic viscosity i/ = ry/p is then: 
I' — 2.7 X 10^^ cm^s"^, tens of orders of magnitude higher 
than in usual fluids. 

As mentioned in sect. El observations of the filaments 
perpendicular to the Galactic plane can be interpreted 
as tracers of a coherent, vertical magnetic field. Magnetic 



^ The viscosity of a hydrogen plasma would be 8 times 
higher. 



fields make the stress tensor anisotropic: the perpendicu- 
lar mean free path is reduced to the Larmor radius, so that 
the diffusive properties are very different for the parallel 
and for the perpendicular direction. When applied to the 
transport of perpendicular momentum, it is found that the 
shear viscosity implying terms in diVj is reduced by the 
magnetic field whereas the bulk viscosity im plying terms 
in diVi remains unchanged l)Braginskiilll965|) . In a general 
manner, five different coefficients determine the stress ten- 
sor of a magnetized plasma. All of them can be expressed 
from 770 as negative powers of flcT where ilc is the cy- 
clotron frequency and r is the collision time. When the 
field is strong enough, all of them vanish b ut rin. The col- 
lision time of the 8 keV helium plasma is ( Belmont et alJ 



T 46 X lO"" yr 



UkeVy V 



.035 cm" 



(3) 



Even for the lowest estimates of the magnetic field 
strength {B ^ l^G), the ratio of the highest coefficient to 
?7o is: rji/rio ~ 10~20. The shear viscosity in the Galactic 
center region is fully inhibited by the magnetic field and 
only the bulk viscosity can diss ipate energy. The bulk vis- 
cous stress is (lBraginskiilll965|) : 



^0 



VD-d\\D e 



(4) 



where the subscript || notes vector components along the 
field, e|| is the unit vector pointing along the field direc- 
tion, and 



D^V.v- 3d, 



(5) 

is strongly related to the fluid compression. The associated 
power locally dissipated by the bulk viscosity is simply: 



(6) 



The implications of this viscosity are quite differ- 
ent from those of the shear viscosity. In particular, the 
Reynolds number may not be pertinent to characterize 
the viscous regime of a flow. For the typical cloud velocity 
Vc and radius r^, the Reynolds number would be: 



7^e 



0.06 



5 pc 



100 km s~^ 
P 

2. 1025 g cm- 



3 keV / ^ ^ 



which would indicate an extremely viscous plasma. 
However, this number was defined as the ratio of the iner- 
tial force over the viscous stress for a shear viscosity. The 
latter depends on terms in diVj which are of the order 
of Vc/vc whereas the bulk viscosity implies terms as diVi 
which can be much smaller. As a consequence, the dimen- 
sionless number corresponding to the bulk viscosity, TZbi 
is larger than the Reynolds number. For the cloud motion 
in the Galactic center, we can only ensure that TZb > -06 
and as a result, the fiow around molecular clouds must 
be studied in details to be able to determine its viscous 
regime. 
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3.2. First example 

To better characterize the way the bulk viscosity acts on 
the flow, we perform here a 2-D analysis. We study the 
flow around an infinite cylinder whose axis is parallel with 
the ambient magnetic field. For a 2-D fiow invariant along 
the field direction, d\\v\\ = and D reduces exactly to 
the flow divergence: D — 'V±.v±. A solution for the flow 
is presented as online material in appendix ^ and it is 
found that, in this geometry, the magnetic field only acts 
as an additional pressure. The solution is thus the same 
as for a neutral gas, except that the wave propagation 
speed is the transverse fast magnetosonic velocity vp- In 
that sense, the flow can be regarded as only resulting from 
standing and propagating fast perturbations. In the low- 
P hmit {B > 100/zG), the dispersion equation for fast 
waves is: = k'^v\. The cloud only excites waves with 
uj ~ Vc/rc and k ~ I/tc 
hot and faint phase is: 



The Alfven Mach number in the 



niA = 0.17 



10-25gcm 



-3 



1/2 



B 



.ImG 



lOOkms" 



so that the clouds of the Galactic center region are very 
sub-Alfvenic. The cloud motion thus does not efficiently 
excite fast waves and we find that the flow around the 
cylinder is mostly incompressible. The flrst compressible 
contribution only appears to the order {vc/v-p)'^: 

V.v - vl/r^Vc (8) 

This leads to define a 2-D bulk Reynolds number: 

7^2z,=3^f^)" (9) 

We see that the stronger the magnetic field, the higher 
the bulk Reynolds number. It has already been noted that 
the magnetic field inhibits the shear viscosity; we see here, 
for this 2-D case, that it also lowers the efiiciency of the 
bulk viscosity. Typically, the bulk viscosity is reduced by 
2 orders of magnitude. For the properties of the Galactic 
center gas, the 2-D bulk Reynolds number is: 



5 pc/ VlOO km s 



B 



.1 mG 



8 keV 



-5/2 



(10) 



The 3-D situation is much more complex and it is 
not obvious how the 2-D conclusions can be adapted. 
First, in a fully three dimensional problem, the viscos- 
ity stress depends on the parallel motion of the fiow. If 
<< 'V±.v±, the 2-D analysis holds and the above di- 
mensionless number can be used to determine the viscous 
regime. On the contrary, if ~ v^/rc, the real Reynolds 
number corresponds to the bulk Reynolds number. And in 
between, when d\\v\\ ~ V^.t;^, the bulk Reynolds number 
is intermediate. Thus the regime can be highly or weakly 
viscous depending on the exact flow properties. Second, 



whereas in this 2-D example only the fast waves seems to 
play a role, other modes, namely the Alfven and slow ones 
make the flow properties in 3-D very different from this 
simple example. As a consequence, a precise description 
of the 3-D flow around a spherical cloud is required to 
estimate the power dissipated in the plasma. 

4. Friction 

In this section, we investigate how the viscosity can dis- 
sipate the energy of the molecular clouds in the plasma. 
For the sake of simplicity, we assume that the viscosity is 
not strong enough to modify the flow. We can first char- 
acterise the flow around one single cloud for an inviscid 
plasma. The dissipation is then estimated by using this 
inviscid solution and the viscosity properties. 

4.1. The clouds wake 

In this section, we refer to already existing literature to 
characterise the wake of molecular clouds flowing in the 
Galactic center region. 

It is often assumed that strong currents can exist at 
least at the cloud surface. If the non-thermal filaments 
result from the field-cloud interaction, the clouds must 
be conducting, at least at their surface. Several processes 
can ionize the clouds enough to achieve a high conduc- 
tion, such as the irradiation by local sources (stellar cl us- 
ters for example, iMorris fc Serabv'nlll996t iMorrij 19961) or 
the Alfven critical velocity ionization effect taale ev et all 
11984 iMorris fc Yusef-Zadeh 1989) . In the following, we 
will assume that they are perfect conductors. We are thus 
interested in the wake of a conductor moving in a magne- 
tized plasma. 

The behavior of electrically conducting bodies embed- 
ded in a flowing plasma has already be en investiga t ed but 
in different astrophysical conditions. iDrell et alJ l|l96,■^^ 
were among the first to address this question. Their goal 
was to explain the anomalous drag experienced by the 
Echo I artificial satellite in the earth magnetosphere. They 
introduced the important concept of 'Alfven wings', which 
are standing Alfven waves attached to the conductor. 
Many studies were also dedicated to the Alfvenic wake of 
the Jovian satellite lo, which is thought to have a goo d 
electrical conductivity ijGoldreich fc Lvnden -Bcl? 'l9 69|). 
An exact solution was found for the wing s (^Neubauen 
ITgillSouthwood et a].ll98nHNenbaneTll99l). which gives 
predictions for fiow and field perturbations in excellent 
agreement with observations from Voyager 1 ( Acuna et alJ 
Il98lt iBelcher et al.lll98ll iBarnett fc OlbertI 11984 - Also, 
the possible contribution to the wak e from non Alfvenic 
perturbations was early suggested llChu fc GrossI Il96fil: 
IWolf-G]a,drow et alJlTgH iT.inker et a,l.lll988 1 . B^consid- 
ering conservation laws, IWright fc Schwarta l|1990l) have 
shown that indeed the Alfven wings cannot exist on their 
own: other disturbances must exist. In the MHD approxi- 
mation three different modes can propagate in the plasma: 
the Alfven mode but also the slow and fast magnetosonic 
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Fig. 1. Schematic view of the MHD wake a a spherical 
cloud. 



ones. The motion of a moving body excites a priori the 
three modes so that the wake can be interpreted as con- 
stituted by three components. When the only source of 
disturbance is assumed to be the body itself, all the waves 
excited in its vicinity propagate away. As their propaga- 
tion properties are different (see below) , they can spatially 
separate and the three components can be easily identi- 
fied at large distance. The general picture of such a wake is 
presented in fig. ^ Numerical simulations present results 
very similar to this simple sketch. They confirm in par- 
ticular the main properties of the Alfven wings and show 
evidence for the slow and fast magnetosonic perturbations 
in the wake (Linker ct al. 1988, 1991, 1998). However, be- 
cause of their higher complexity compared to the Alfvcnic 
ones, the non Alfvenic disturbances were very little stud- 
ied. 

Since the three MHD modes have different properties, 
they are not subject to the bulk viscosity with the same 
efficiency (see sect. 14 .21) . It is thus important to estimate 
their respective amplitude. The coefficients which deter- 
mine the relative contribution of the three modes to the 
wake depend on the exact conditions at the cloud surface, 
such as the cloud speed, its shape, its conductivity, the 
magnetic flux its carries etc... The characteristics of the 
cold molecular clouds in the Galactic center region are 
different from those of lo and the artificial satellites and 
they may not excite the three modes in the same manner 
as satellites do. 

The relative velocities in the Galactic center are in par- 
ticular very different from those of the satellites wake pre- 
viously studied in the literature. The motion of satellites 
like lo or Echo I is rather fast compared to the wave speed: 
t'c Cs ~ WA, whereas we are interested in very subsonic 
motions: Wc < Cg ^ ua- For this reason, our conclusions 
cannot be compared directly to the results of numerical 
simulations that exhibits solutions for low Mach numbers: 
rriK = .5 and = .5- 1 (Linker et al. 1988.. ,1991, 1998). 

In the followin g, we will extens ively refer to the so- 
lution proposed bv iNe^ibaiied 1 1980 ^. which describes the 
wake of a magnetized satellite carrying the same magnetic 



field as in the surrounding medium. The magnetization of 
molecular clouds is not well understood. Their conduc- 
tivity can be expected to be high enough to significantly 
prevent the magnetic field to diffuse in or out of the cloud 
on the crossing time of the cloud so that the magnetic 
flux threading the cloud remains constant on this typical 
time scale. The present magnetization thus depends on the 
poorly constrained initial magnetization and on the cloud 
history on long time scales. It is often argued that the field 
inside the clouds is made of external field lines which are 
tangled by the turbulent velocity of the inner clumps. The 
clouds magnetization must thus be a finite fraction of the 
external one. On the one hand, the cloud ionisation may 
be only partial so that the diffusion time may be shorter 
than the cloud life time. In such case, the magnetic field 
lines may have diffused on long time scales as the cloud 
moves, so that internal and external fields could now be 
comparable. On the other hand, it has been proposed that 
the external magnetic pressure could be responsible for 
the clouds confinement, pointing to a weaker cloud mag- 
netization. The case where the cloud does not carry any 
field is very complex and has not been studied in details. 
Nevertheless, it has been proposed that this case would 
lead to an exact solution different from the wake with a 
strong magnetiza tion, but with similar general properties 
l|Neubaueii Il998l sect. 5.1.2). For the sake of simplicity, 
we will thus assume that the clouds carry fields lines with 
comparable strength to those of the external medium. 

In the next three subsections, we describe the proper- 
ties of the different modes and their respective amplitude 
for the subsonic motion of a conducting body in a strongly 
magnetised plasma (/3 << 1). We stress in particular the 
compression and the parallel velocity which govern the 
dissipation efficiency. From now on, we refer to the ver- 
tical direction as the field aligned one and the horizontal 
plane as the plane perpendicular to B and including the 
cloud-field relative velocity. The subscript || refers to the 
direction parallel to magnetic field, whereas the subscript 
_L refers to the components perpendicular to the mean 
field. We use the frame moving with the cloud. 

4.1.1. The Alfven wings 

The Alfvenic component of the wake has been most stud- 
ied. Since the MHD equation admit a non-linear solution 
for the propagation of Alfven waves , an exact non-li near 
solution has even been presented bv lNeubaueil 1 198r^ . 

The Alfven waves are excited at the cloud surface by 
the transverse velocity and the conditions on the electric 
field; namely, the electric field tangent to the cloud sur- 
face has to vanish in the comoving frame. Their prop- 
agation is strictly along the field lines, which act as a 
guide. Thus their amplitude does not decay by geometric 
effects, as it would do if they were emitted isotropically. 
As the Alfvenic perturbations travel vertically, the flow 
drags them horizontally. As a result, the Alfven waves ex- 
cited in the cloud vicinity are always located along the 
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Alfven characteristic, forming the so called 'Alfven wing'. 
There are two of them, one on each side of the cloud, cor- 
responding to the two possible directions along the field. 
Each wing axis corresponds to an invariant direction. They 
are inclined at an angle q;a — tan^^rriA from the verti- 
cal direction. The clouds are very sub-alfvenic, so that 
the Alfen wings are almost vertical: ofA ~ -02 — 0.2 for 
B w .1-1 mG. 

By means of the Alfven waves, the condition E' = 
at the conductor surface propagates along the field lines, 
so that the comoving electric field vanishes in the whole 
flux tube threading the cloud, forming a cylindrical wing. 
Inside the wing, the magnetic field is aligned along its axis 
and the perpendicular velocity vanishes. The cloud thus 
entrains the whole flux tube and its full content along with 
it. As mentioned earlier, we assume in first approximation 
that the field threading the cloud is comparable with the 
external field. The wing size is thus of the order of the 
cloud diameter. Outside the wing, the stream and exter- 
nal field lines avoid the cylinder: this corresponds to per- 
turbations decaying as 1 jr^ where r is the distance to the 
wing axis. The surface of the cylinder is characterized by 
a current sheet which closes the currents genera ted within 
the co nducting cloud. In the solution given bv iNeubauerl 
l)l98fl) . each wing corresponds to a Poynting flux of: 

i^A = T^rlpvlvp, (11) 

As the stream avoids the whole flux tube, the per- 
turbed velocity is comparable with the cloud velocity: 
Sv±/vc ~ 1 inside and in the vicinity of the cylinder. This 
wing is fully incompressible: 5p/p = 0, even at the non- 
linear stage. The linear solution has no parallel velocity, 
so that the linear Alfven wing has D — and no dis- 
sipation. However, the non-linear solution allows a small 
parallel velocity which allows a finite D and thus a small 
dissipation. When expanded to the first non-linear order, 
it is found that: 

Ull/Wc^TlA (12) 

As we will see, the Alfven wing derived by Neubauer 
is the main contribution to the wake, but it does not fully 
satisfy the boundary conditions at the cloud surface, so 
that contributions from the other modes are needed. 

4.1.2. The slow wings 

For instance, the parallel velocity found in the non-linear 
solution for the Alfven wings induces a net outflow. As 
the cloud cannot provide this material, this mass flux 
must be balanced by other perturbations. The respective 
role of the fast and slow modes is not well understood. 
IWright fc Schwartz! lll99(l) suggested that the matter sup- 
ply can be at least partially attributed to the slow pertur- 
bation. 

The general propagation of slow magnetosonic modes 
is more complex than that of Alfven ones but, in the low-/3 
limit, the slow waves are also guided along the field lines. 



As Alven perturbations do, they are thus able to form 
wings on both sides of the clouds. The slow modes how- 
ever propagate with a velocity close to the sound speed, 
slower than the Alfven one, so that the slow wing an- 
gle tts — tan^^ TOg is larger. For a typical sound speed of 
1000 km s^^, as = 0.1. For mG fields, the slow wing sepa- 
rates from the Alfven wing at an altitude h « 2rc/ms, i.e. 
^ 100 pc. For lower fields, they separate even farther, but 
when the field is in equipartition or weaker, slow waves do 
not strictly propagate along the field lines, so that they 
do not form a truly infinite wing. As long as the field is 
not too weak, the global shape is however similar. Since 
we consider a 70 pc height region, the two wings are thus 
superimposed in the whole emission region for any field 
strength we consider in this paper. In the linear approach, 
we can however study them independently. In the follow- 
ing, we will assume that the parallel velocity of the slow 
wave balances the parallel velocity of the Alfven wings. 
The polarization relation of slow waves gives: 

= -fi^n (13) 

The slow wings are almost vertical so that k\^lk^ « 
TOg << 1. As a result, the transverse perturbed veloc- 
ity remains very small in comparison to the relative ve- 
locity and to the Alfvenic transverse velocity. From the 
slow mode properties, we also get the perturbed density: 
bpj p ~ mpjrig,. Since the motion of molecular clouds is 
subsonic and sub- Alfvenic, the compression is very weak: 
Spj p « 1. In the low-/3 limit where the slow waves form 
a wing, the quantity D in the slow wings is dominated by 
the contribution of the parallel velocity: D « —2k\\v\^. 

4.1.3. The fast perturbation 

It is likely that fast perturbations are also excited. It has 
already been menti oned that the Alf venic perturbations 
in the solution from INeubauerl l|l98fll) extend beyond the 
cylinder frozen in the cloud. Nevertheless, the Alfvenic 
perturbations directly excited by the cloud surface can 
only propagate in the flux tube threading the cloud and 
Alfvenic perturbations propagating on external field lines 
cannot be excited by the cloud itself. They must be non- 
linearly excited by an other mode which would directly 
be generated by the cloud. The slow mode is guided by 
the field lines and can not contribute to this excitation, so 
that only a fast perturbation extending around the cloud 
can be responsible for the external structure of the Alfven 
wing. This shows that the linear analysis must be handled 
with care, but also that there must be a fast component 
in the wake. 

The general propagation of fast magnetosonic modes 
is complicated too but, in the low-/3 limit, the fast waves 
propagate isotropically, as simple sound waves at the 
Alfven speed. They thus cannot form a wing and their 
amplitude decays with the distance to the source. As a 
result, the fast perturbations excited by the cloud remain 
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localized in its vicinity, just as in the classical hydrody- 
namical wake. 

Because of the spherical propagation of fast waves, 
it is not easy to derive an estimate for D in a fully 3-D 
case. We can however try to rely on the 2-D results. We 
have found that, for a subsonic and sub-Alfvenic motion, 
the cloud cannot excite waves that satisfy the dispersion 
equation of fast modes, so that the compression is very 
weak: Sp/p ^ m\. We moreover note that this is also the 
case for a general 3-D wake in a non magnetised gas. The 
same conclusions must apply in 3-D for our MHD flow. 
We will thus use this relation for our following estimates. 
In the low-/3 limit, the polarization relations of the fast 
plane waves give: fc||W|| = (3k±v±. Assuming it is still the 
case for spherical waves, we find that the parallel velocity 
has only a very small contribution in the quantity D, so 
that it is dominated by the perpendicular flow divergence: 
D ^ 'V±.v±. If not, this result can be used as a minimal 
estimate of the dissipation in the fast perturbation. 



4.2. Viscosity efficiency 

Because of their different properties, viscosity acts differ- 
ently on the 3 components. In the absence of any other 
cooling mechanism, the dissipation efficiency can be de- 
termined by comparing the power dissipated by viscosity 
and the X-ray emission in the whole central region: 

" ^ (14) 



L w 4 X 10^"^ erg s"^ 



70 pc 



Here we specify this efficiency for the three different wake 
components. 

4.2.1. The Alfven wings 

As mentioned earlier, the Alfven waves are incompressible, 
and have no parallel motion in their linear form. Hence, 
in the Alfven wing, D = and there is no dissipation due 
to the bulk viscosity. To lowest order, the energy in the 
Alfven wings thus leaves the central region without being 
damped nor heating the plasma. However, we find that the 
integrated Poynting flux for both wings of the clouds 
in the Galactic center is typically one order of magnitude 
higher than the X-ray luminosity: 



Fa = 2. X 10^^ erg s 



VlOO 



p 



5 pc 

1/2 



B 



,1 mG 



(15) 



2.10-25 g cm-3 



, 100 km s" 

It is thus natural to wonder if secondary effect could dissi- 
pate a small fraction of this large energy. For 100 pG fields, 
10% would be required whereas for mG fields, only 1% is 
needed. In general, the damping of Alfven waves is not 
very efficient, but different processes can dissipate a small 
fraction, as for example small magnetic irregularities, den- 
sity gradients... Here, we concentrate on two possible vis- 
cous mechanisms: the damping of the non-linear parallel 



velocity and that resulting from a small compressibility of 
Alfven waves propagating along magnetic field lines which 
have a large scale curvature. 

As mentioned previously, lNeiibaiie'3 l|l980l) could ex- 
hibit an exact non-linear solution for the Alfven wings. 
This solution, although it is incompressible, involves a 
small parallel velocity which allows dissipation by the bulk 
viscosity. From eq.^jwe derive D 3m\vc/rc. As long as 
the viscosity does not dissipate the whole energy flux, we 
can estimate the energy dissipated along the wing up to 
an altitude of 70 pc by summing the local dissipation rate 
over this height. For 100 clouds, we get the dissipation: 



Qai = 2. X 10^^ erg s"^ 



100 



70 pc 



2.10-25gcm-3 



100km s" 



B kcV 

B 
.ImG 



5/2 

(16) 



Since the non linear parallel velocity is proportional to 
Vcl'^A^ the magnetic field strongly inhibits this non linear 
dissipation and for fields above equipartition, this power is 
less than the radiative cooling. However, we find that for 
fields lower or equal to 100 /iG, this energy is comparable 
with the radiative losses and can therefore account for the 
plasma heating. 

On the other hand, the observed non-thermal filaments 
seem to be slightly bent, which may trace a large scale 
curvature. From observations, we estimate the radius of 
curvature being of the order of 100 pc. This is consistent 
with the hypothesis that the large-scale vertical magnetic 
field might be due to currents in the molecular torus at 
150 pc. In such conditions, the usual description of plane 
waves can not be applied exactly. In the limit where the 
radius of curvature is large in comparison with the wave- 
length, modes can be identified with the usual Alfven, 
slow and fast ones, but with slightly different properties. 
For instance, Alfven waves become weakly compressional: 
any Alfvenic displacement perpendicular to a curved mean 
field experiences a compression or a depression depending 
whether it points to or away from the curvature center. 
This is illustrated on fig. El This property can be inter- 
preted as a coupling between the slow and Alfven waves 
in a curved magnetic field (S outhw ood & Saunders 198^. 
From fig. 13 we can estimate the divergence of any Alfvenic 
perturbation to be of the order of: 



(17) 



where Vj. is the perturbed velocity component projected 
along the curvature direction and is the local radius of 
curvature. Since they are compressible, Alfven waves in a 
curved magnetic field can be damped by the bulk viscos- 
ity, even at their linear stage. A more detailed treatment is 
presented in appendix (online material) where we show 
that the quantity D responsible for the dissipation is ac- 
tually: 



Da2 



2vc 
Rc 



1 



(18) 
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4.2.2. The slow wings 

We have assumed the parallel velocity of the slow wing to 
be similar to that in the Alfven wings. By the polarization 
relations of slow waves, this fully determine the slow wing 
amplitude and gives in particular D ~ niArrisVc/rc- When 
/? < 1, this is larger than in the Alfven wings because 
the slow wings arc more inclined. To estimate the power 
dissipated in the whole central region, we can again sum 
the local dissipation rate over a height of 70 pc, and we 
find that the 100 clouds of the Galactic center can provide 



Fig. 2. Illustration of the compression associated with the 
propagation of Alfven-like perturbation in a curved mag- 
netic field. The field lines are represented by the parallel 
curves and define a local Alfven velocity ua, -Rc is the 
curvature radius and Vr is a transverse perturbation along 
the curvature direction. For straight lines, the Alfvenic 
transverse perturbation has no compression, whereas in 
a curved geometry, the volume of small elements of flux 
tubes changes with the distance to the curvature center: 
any displacement towards the curvature center drives a 
compression and any displacement away from the curva- 
ture center drives a depression. 



In the low-/? limit, the corresponding dissipation for Nc 
clouds is: 



Qa2 = 3. X 10^^ erg s'^ ( ^ 



100km s"V V5pc 



(ML) 



5/2 



70pcy \8keV J 



(19) 



Vioopcy 



1-3- 



Wc sec that when v\/c^ ~ 3, i.e. when B w 300 iiG 
there is no dissipation. This is however a limit between to 
domains where the dissipation is significant. On the one 
hand, for stronger fields, the dissipation tends to a finite 
value which depends on the curvature radius. We find that 
a curvature radius Rc ^ 100 pc is required for the damp- 
ing to balance the radiative cooling in a mG field, which 
may be consistent with the observations. And on the other 
hand, for weaker fields, the dissipation is quickly compa- 
rable with the radiative losses and can become orders of 
magnitude stronger if the field strength is smaller than 100 
/xG. Obviously, the bulk viscosity cannot dissipate more 
that the Poynting fiux. When our estimate exceeds the en- 
ergy fiux of the Alfven wings, the viscosity starts changing 
the flow properties and our estimate fails. A full treatment 
of the viscosity is then required. However, it is likely to 
give a saturated solution where the Alfven flux is fully 
dissipated and this energy is more than needed to balance 
the radiative cooling. 



Qs = 3 X 10^'=' erg s 



Nr. 



VlOOy \70 pcj \8 kcV 
p \ I B 



5/2 



10-25 g cm-3 

6 



.1 mG 



100kms"V V 1000 km s"^ 



(20) 



We find that the power dissipated in the slow wing is 

one order of magnitude lower than the radiative losses. To 
reach a balance, weaker fields or a higher accretion rate 
are thus needed. 



4.2.3. The fast perturbation 

Fast waves arc compressible. They coiild then be expected 
to be responsible for a strong dissipation. However, as we 
have seen, since the motion is subsonic, the compression 
remains weak whereas the velocity perturbation is large. If 
Spii < V_L.t»_L, the results found in the two-dimensional 
analysis are still valid in 3 dimensions and D ~ m\vc/rc- 
This value is comparable with the non-linear Alvcnic D. 
However, contrary to the latter, the fast perturbation is 
located in the cloud vicinity. The total efficiency is thus 
lower. We can estimate the total power dissipated by mul- 
tiplying the 2-D results by the height of the cloud 2rc. 
This leads to: 



/2 



Vr \ I B 



10-25gcm-3/ V 100km s"V V-lmG 



(21) 



This is two orders of magnitude too small to contribute 
significantly to the dissipation. As we have noted earlier, 
if the parallel velocity plays a more important role in the 
3-D situation, then this might underestimate the dissipa- 
tion. Moreover, even if the fast perturbation amplitude 
and the associated dissipation decay with distance to the 
cloud, the dissipation integrated over far spherical shells 
might remain significant. This would lead to using a larger 
scale height. It is thus possible that the total dissipation in 
the fast perturbation is stronger than the above estimate. 
A more detailed treatment is therefore required but the 
dissipation is unlikely to dominate over the Alfvenic dis- 
sipation, so that we do not discuss this possibility further 
on. 
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Fig. 3. Total dissipation versus magnetic field strength for 
different radii of curvature Re- The power is normalized to 
the X-ray luminosity in the Galactic center region, so that 
1 is the power required to balance the radiative cooling. 
Contributions Qai, Qs and Qy for the power dissipated 
in the Alfven wing by non-linear effect, in the slow wing, 
and in the fast perturbation respectively, are also plotted, 
although the estimates for the slow and fast perturbations 
lose accuracy below 100 ^G. 

5. Discussion and conclusion 

As a conclusion, we find that the dissipation is different in 
each component of the wake and that the global efficiency 
depends on the field strength. For all strengths, dissipa- 
tion within the fast and slow wings is negligible and most 
of the dissipation occurs in the Alfven wing. We find that 
the dissipation resulting from the curvature of field lines 
dominates for any field. Figure O shows the total dissipa- 
tion and contributions for each wing, as a function of the 
field strength. Three domains can be observed: 

— For low fields: B < 100/iG, both the non-linear and 
the curvature damping dissipate enough power to bal- 
ance the X-ray luminosity. The dissipation due to the 
curvature of field lines seems to be dominant, for any 
reasonable radius of curvature. In this strength range 
for the magnetic field, the dissipation in the Alfven 
wings can be responsible for the plasma heating. 

— For intermediate fields: B « 0.2-0.5 mG, the dissipa- 
tion is more than one order of magnitude lower than 
the losses by radiation. 

— For strong fields: B > 1 mG, we find that none of 
the components can dissipate enough energy to bal- 
ance the radiative cooling in a straight geometry. 
Nevertheless, a large scale curvature with Rc < 100 
pc would be responsible for a sufficient dissipation. 

We see that except for the indirect dissipation, the 
magnetic field strongly limits the efficiency, favoring low 
fields. However it must be kept in mind that, for too weak 
fields, some of our results are not valid anymore. Indeed, 
we used the low-/3 plasma limit to allow a very simple 



characterization of the slow and fast perturbations. In this 
limit, the slow modes are guided and the propagation of 
fast sonic modes is isotropic. When B < 100/j,G, the field 
is below equipartition and the properties of the fast and 
slow waves are more complex. However, it is likely that 
as long as the field is not too weak (B > 10 fJ.G), our 
related conclusions remain valid in order of magnitude. 
The description of the Alfven wing does not depend on 
any assumption on the field strength, so that even for 
weak fields, our related conclusions still hold. 

Our model however strongly depends on the field ge- 
ometry, which is debated too. The field has for example 
been suggested to be tur bulent and mostly to roidal as 
in the rest of the Galaxy In such 

case, our results cannot be applied to describe the inter- 
action of the cold molecular clouds with the magnetic field. 
Nevertheless, the viscosity may still play a role. In a verti- 
cal field, most of the energy can fiow away as Alfven waves 
without being dissipated. In a turbulent medium, it has 
to stay within the plane. Coupling with other modes fi- 
nally generates compressible modes which then have to 
be damped by viscosity. For these reasons, it is likely 
that, even in a turbulent field, the motion of the molec- 
ular clouds can participate to the plasma heating. The 
dissipation in a turbulent medium requires a dedicated 
investigation and it is not excluded that it might be even 
stronger than in the straight geometry. 

These results also depend on the statistical properties 
of the clouds, which are poorly constrained. In particu- 
lar, we have assumed clouds of spherical shape, which is 
obviously not accurate for the clouds in the Galactic cen- 
ter. The dissipation generated by their extended shapes 
and fractal structure cannot be easily derived, but local 
enhancements of the fluid compression are expected so 
that the global heating could be stronger. Following the 
2-D analysis presented in appendix we computed nu- 
merically solutions for non-cylindrical clouds flowing in a 
viscous medium. We found that multi-pole shaped clouds 
may dissipate several times more than a simple spherical 
cloud. The real structure of the clouds is far more complex 
than multi-poles and could result in an even more efficient 
dissipation. 

The dissipation rate in the non-linear Alfven wings, the 
slow wings and the fast perturbation also strongly depends 
on the cloud velocity: Q oc v^. Any uncertainty on the es- 
timate of this parameter can therefore have a strong effect 
on the order of the dissipation. In particular, in our esti- 
mates, we used an average value for the velocity whereas 
we should integrate the power resulting from each cloud 
over the velocity distribution. Such an integration would 
favor the fastest clouds and give a higher total dissipation 
rate. As the cloud velocity dispersion is large, this effect 
may be important enough to make the dissipation in the 
slow wing competitive with the radiative losses, even for 
intermediate field strengths. 

We also note that the accretion of clouds from the 
molecular ring at 180 pc is not thought to be continu- 
ous. Some clues seem to indicate that it could be strongly 
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intermittent with catastrophic events every ~ lO'' yr 
llStark et al J 12004 . The number of clouds and their ve- 
locity are expected to be higher during these events, pro- 
viding a stronger heating. Because the radiative cooling 
time is very long (> 10^ yr) compared to the inter-event 
time, such a past heating would not be distinguished from 
a present one. 

A very interesting property of this heating mechanism 
is that the viscosity regime seems to correspond exactly 
to the maximal efficiency for dissipation. We have pre- 
sented here a MHD analysis of the cloud-field interaction. 
This holds as far as the collision time is shorter than the 
typical time scale of the problem. From eq. El it is seen 
that when the friction heats the plasma, the latter be- 
comes more and more viscous. But meanwhile, the col- 
lision time increases (eq. The plasma thus reaches a 
temperature where the regime becomes coUisionless. If the 
temperature keeps gro wing, the c o llision s become so rare 
that the results from iBraginskiil 1(19614 ) fail. The wave- 
particle interaction must then be studied in a kinetic for- 
malism, leading to damping by magnetic pumping; but 
as the viscosity basically results from collisions, it must 
drop in the coUisionless regime. The maximal dissipation 
is thus reached when the collision time is of the same or- 
der of magnitude as the typical time scale of the problem. 
In the case of the cloud motion, this time is the cloud 
crossing time: At = Tc/vc « 48 x 10'^ yr, showing that 
the conditions are optimal for a maximal dissipation. This 
naturally suggests a self-regulation mechanism where the 
viscous friction heats the plasma to the temperature that 
makes the regime weakly collisional. Then, the efficiency 
drops and the temperature saturates. Given the strong dc- 
pendance of the dynamical viscosity and the collision time 
with temperature, this gives a precise saturation temper- 
ature, which could be ksT—S keV. 

The overall interaction of the moving cloud with the 
ambient field exerts a drag on the cloud and must slow 
it down, causing it to accrete toward the Galactic cen- 
ter. In the processes we have presented here, most of the 
energy taken from the cloud goes away as Alfven waves, 
of which only a fraction is dissipated by the viscous pro- 
cesses we have described. The rest of the energy leaves the 
X-ray emitting region without contributing to the heat- 
ing but does participate to the drag on the cloud. This 
phenomeno n successfully exp lained the unexpected drag 
on Echo I l|Drell et al.lll965() and may be significant for 
the molecular clouds too. We estimate this energy to be 
2 ~ 10'^^ erg s""'^ per cloud. This gives a typical drag- 
induced radial velocity of Vr ^ 1 km s~^ (B/l niG) and 
an accretion time of ~ 10^ yr {B/1 mG) . For strong 
fields, these results are comparable with the other accre- 
tion mechan isms such as the dynamical friction with stars 
l|Stark et al.Ll991) . They are also consistent with the ac- 
cretion rate needed to sustain a constant accretion from 
the external mo lecular torus at 180 pc and explain the 
mass profile ([Morris fc Serabvnlll99(it) . The pervasive field 



in the Galactic center may thus significantly contribute to 
driving the gas accretion. 

As a conclusion, we find that in spite of the strong 
vertical field, the viscosity in the Galactic center region 
can dissipate enough of the kinetic energy of the molec- 
ular clouds that orbit in the CMZ to balance the X-ray 
luminosity. This process not only provides the necessary 
power, but is also able to bring the gas to the observed 
temperature of 8 keV, which turns out to be the one where 
this process is most efficient. It should also contribute to 
the drag causing the clouds to accrete at this intermediate 
scale between the large-scale, bar-dominated fiow and the 
smaller scale of the central black hole vicinity. 
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Appendix A: The wake of an infinite cylinder 

In this appendix, we study the two-dimensional subsonic, 
sub-Alfvenic MHD flow around a cylinder for any viscous 
regime. This allows us to proceed by expansion starting 
from an incompressible solution. 

A.l. Equations 

We study the flow around a vertical cylinder of radius Tc. 
The magnetic field far from the cylinder is assumed to be 
parallel with the cylinder axis. The flow at inflnity has a 
constant velocity Vc and is perpendicular to the axis. 

In this simple two-dimensional case, the viscous stress 
and the Lorentz force write very easily. On the one hand, 
the viscous stress is proportional to the compression of the 
fluid: 

F^ = ^V±{r]oV±.v±) (A.l) 

On the other hand, in this geometry, there can be no per- 
pendicular perturbed magnetic field so that the field lines 
remain vertical. They only compress and avoid the cylin- 
der. As a result, the magnetic field only acts on the plasma 
dynamics by its pressure that is proportional to the den- 
sity. 

Finally, the system constituted by the mass conserva- 
tion and the equation of motion reads: 

V.{pv) = 

p{v.V)v = -vl,Vp+^V{r]V.v) 

where = Cg -|- v\ is the fast magnetosonic velocity. It 
is non linear and cannot be solved at once but further 
assumptions can be made. As noticed in the main body of 
the paper, the motion of the molecular clouds is subsonic 
and sub-Alfvenic. We can thus restrict our study to this 
case and order the equation with the small parameter: 



where is the fast magnetosonic at infinity upstream. 
Prom now, we use dimensionless quantities. Velocities, 
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density and viscosity are normalized by their value at in- 
finity upstream whereas distances and gradients are nor- 
malized by the cylinder radius. Then, the equation of mo- 
tion reads: 



ep{v.V)v ^ -vlyp 



1 



V (77V. v) 



(A.3) 



where we have defined the bulk Reynolds number: 

Tl2D - ^ (A.4) 

We can now order these equation with respect to e. 

A. 2. Lowest order: incompressible flow 

To lowest order, the incrtial terms do not appear. By in- 
tegrating once the stream aligned component of the Eulcr 
equation, we get: 



ho{pa) 



Vo 



(t;o-V)po = cste 



(A.5) 



where /ig = ^i^^ j -I- 2f^ is an extended enthalpy for the 
fluid. Equation IA.5I is a first order differential equation 
for po and there is no solution that does not diverge but 
the incompressible one: po = 1; V.^o = 0. In cylindrical 
coordinates (r, 6) and with regular boundary conditions 
{vr = at the cloud surface and Voo — Vc) we find: 



cos 9 



' sm f 



1 



(A.6) 
(A.7) 



Figure lA.ll shows this solution where stream lines avoid 
the cylinder without compressing. This solution is purely 
incompressible, so that the bulk viscosity can not dissipate 
energy. 

A.3. First order: compressible solution 

To following order, there are compressible terms. As be- 
fore, we can integrate once the equation of motion to: 



(A. 



Contrary to the zeroth order equation, the source vanishes 
far from the cylinder. We can thus find a non diverging 
solution. The perturbed density is then found by integrat- 
ing this equation over the stream lines. In the limit cases 
where 7?.2_d << 1 or TZ2d >> 1, the solution can easily be 
found whereas it is easier to make the integration numeri- 
cally otherwise. Figure TKA\ gives the perturbed density pi 
for TZ2D = 1- This solution is compressible so that, to this 
order, the bulk viscosity dissipates energy. The dissipation 
is proportional to the square of the fiow divergence which 
is easily derived from the perturbed density by the mass 
equation: 




Fig. A.l. Total dissipation per unit height in the MHD 
flow around an infinite cylinder for TZ2d- Flow comes 
from the left hand side. The upper, mid and lower pan- 
els show the main velocity, the perturbed density and the 
perturbed flow divergence, respectively. Dashed and plain 
contours are negative and positive contours, respectively. 



Figure TKa\ shows the contours for the flow divergence. 

To get the total dissipation per unit height, the local 
one must be integrated over the horizontal plane {{x,y). 
In a general manner, it depends in the viscous regime and 
can be written 



Q/Qo 



where 



1 



V.vi 



dxdy 



Qo = Parc-§- 

V 



(A.IO) 



(A.ll) 



is the typical dissipation and all the quantities in the right 
hand side are dimensionless. For comparison, the dissipa- 
tion by the usual shear viscosity would be: 



i she 



PorcV^ 



(A.12) 



{vo.V)pi = -V.vi 



(A.9) 



If the viscosity coefficient were equal, the bulk viscosity 
would be lower than that by the shear viscosity, by a factor 
e. Figure IX!2l shows the dissipation by the bulk viscosity 
as a function of the bulk Reynolds number. 

A.4. Conclusion 

We find that the fiow in the subsonic limit can be de- 
composed as a main incompressible solution and a small 
compressible perturbation. The dissipation rate depends 
on the viscous regime. The maximal dissipation occurs 
for the intermediate regime when TZ2D ^ 1- The corre- 
sponding estimate is e = (fc/fm)^ times smaller than the 
dissipation which would occur via shear viscosity without 
magnetic field. 




Fig. A. 2. Dissipation rate as a function of the 2-D bulk 
Reynolds number. Here, Qq = poVcV^/v'^. Points are the 
result of numerical integration of equation IA.8I The two 
asymptotes Q/Qq = 1.837?.2_d and Q/Qn — BAbTZ^^ re- 
sult from its analytical integration in the viscous and non 
viscous regimes. 




Fig. B.l. The equilibrium magnetic field and notation 
definitions 



Appendix B: Propagation of Alfvenic 

perturbations in a curved magnetic field 

In this appendix, we derive the characteristics of the 
Alfven waves propagating in a curved magnetic field. As 
soon as the medium is inhomogeneous, the plane waves are 
not the eigensolutions and the true solutions have often a 
complex behavior. In the limit where the irregularities are 
weak, the solutions are similar to plane waves, but with 
slightly different properties. In particular, we can expect 
the Alfven waves to gain a compressional component, and 
thus to become subject to viscous damping. Here we write 
the full waves equations in a curved equilibrium. Wc then 
expand the equations by considering the curvature of field 
lines as very weak on the wavelength scale. This expansion 
gives the modified waves properties. 

To model the effect of field curvature, it is more conve- 
nient to use a cylindrical symmetry with invariance along 
the vertical axis z and the azimuthal angle 9. As shown in 
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fig. IB. II the equilibrium magnetic field is assumed to be 
purely azimuthal, so that the corresponding current is in 
the vertical direction: 

Bo = Ba{r)ee Jo = -drirBo)e, (B.l) 

r 

In a static equilibrium, the pressure gradient needs to bal- 
ance the Lorentz force: 

drPo + i-^j - (B.2) 

We note that in this geometry, a current-free equilibrium 
is obtained for Bq oc l/r. Nevertheless, for reasons that 
will be explained later, we will use a specific equilibrium 
that makes = v\/r'^ — Bo/Airpor'^ constant. 

B.l. General perturbed equations 

We shall now perturb the system with small quantities 
Pi, Vi and bi. We perform a Fourier series in 9 and a 
Fourier transform in but, since the radial direction is 
not homogeneous, wc keep the radial dependancy: 

Btot - Bo + b^{r)e-'^^'-^<'-''''-'^ (B.3) 
Vtot = + t;i(r)e-*('^*-'"^-'=^^) (B.4) 
Ptot - po+Pi(r-)e-^('^*-'"''''=^^) (B.5) 

From now, we drop the subscript and denote by ' the 
radial derivative. We express the whole wave system as 
function of the three perturbed velocity components Vr, 
vg, and Vz- With these notations, the mass and induction 
equations read respectively: 

iujp/po = v'r + — ve + ikzVz + — [l + r^] (B.6) 
r r \ poj 

iub/Bo = -—V (B.7) 
r 

f I '<Jr f , {rBo)'\\ 

+ \ Vr + — ve+ikzVz + —[-! + ^^^jjee 

The largest difficulty in this geometry is the radial de- 
pendancy that may not be simply expressed as a Fourier 
coefficient. To deal at best with this, it is easier to ex- 
press the system as a function of the radial velocity Vr- 
After substituting p and b in the vertical and orthoradial 
components of the equation of motion, the azimuthal and 
vertical perturbed velocities can be formulated: 

2 { f 2 2 \ I 

+ (-'-^-A-2fcM)^} (B.8) 
ikz ( f 2 2 2 2 \ / 

+ (B.9) 

with 

V^^^-vli^— + klj + —clvl (^— + fc2 j (B.IO) 
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where we have used the fast magnetosonic speed Wp — 
vl + cl 

To get the dispersion equation, these velocities must 
be substituted in the radial component of the equation 
of motion. It can eventually be written in the following 
manner: 



d_ 

dr 



-P[uj, m, kz,r)—rvr 
r or 



+ Q{Lu,m,kz,r)vr = (B.ll) 



where: 

P = 



PO 2 I 2 



m 
r 



2 



2 9 9 



,2„,2 



„2„.2 



(B.12) 



(B.13) 



This equ ation is t h e ge ne ralized H ai n-Liist 
equation llHain fc LiistI Il958t ICxoedbloedl Il97lt 
iGoedbloed &: PoedtsI I 200^ . In a general manner, P 
vanishes at specific radii. Solving this equation with 
boundary conditions thus leads to deal with singularities. 
This effect is not due to curvature but to the gradients 
of Alfven velocity and of the parallel derivative — 
which are not of interest here. Here, we concentrate 
on the curvature effect. To this purpose, we refer to 
the equilibrium quoted earlier that makes ajj^ = fc^w^ 
constant. In this case, the Alfven condition uj'^ — uj\ = 
can be satisfied for any radius simultaneously and the 
singularity disappears. 

B.2. Modified propagation properties 

When oj^ uj\, P vanishes but Q does not. As a result, 
pure Alfven waves are not solutions of the full dispersion 
eauation lB.llI However, we are only interested in a large 
scale curvature of the field lines. We thus only consider 
perturbations whose typical scale length is shorter than 
the radius of curvature. In this limit, the difference be- 
tween the exact Alfven-like solutions and the plane Alfven 
waves must be small. Hence we now expand the equations 
by assuming that the perturbed system varies faster than 
the background quantities, i.e. on length scales smaller 
than r. As we are only interested in Alfven waves whose 
behavior is mostly determined by their parallel wave- 
length, we only assume that the parallel wavelength is 
shorter than r: 



m » 1 



(B.14) 



The expansion is thus done with the small parameter 
l/rn?. The lowest order describes the plane propagation 
of the Alfven waves and higher orders introduce curvature 
effects. P is a lowest order term: P = P_i ~ and Q 
contains both lowest and first terms: Q = Q-i + Qo with 



(B.15) 



The following order term Qo introduce the curvature ef- 
fects and is responsible for a shift in the Alfven frequency. 
To find this shift, we can write oj — w_i -I- loq, substitute 
in eg. IB.lll and identify order per order. To lowest order, 
we find the propagation relation of plane waves: 



9 IIL n 

^-1 = TT^A 



and to first order, we find: 



r por'^^ \ klpQVr J 



(B.16) 



(B.17) 



We can furthermore assume that the radial gradient for 
the perturbed velocity is also stronger than for the back- 
ground profiles. In a WKB approximation, the radial 
derivative can be written drVr ~ ikrVr and: 



1 (Por 



27V 



r Par 



27 



1 



k^ 

k1 



(B.18) 



The shift in the Alfven frequency depends on the back- 
ground profiles. We will however see in the next subsection 
that this shift is not important for our concern. 

B.3. Modified compression and viscosity 

The general compression of the perturbed flow can be de- 
rived from eq. IB. 81 and IB; 



V.v 



m 



— {uj'^—vi]-{rVry-2 



1 



OJ kjvj^ Vr 



(B.19) 



The divergence involves two kinds of terms. Those in 
v'j, corresponds to the usual compression of sonic modes, 
also existing in a straight configuration. For plane Alfven 
modes, these terms vanish. The other terms, in Vr/r, result 
from the field curvature. By substituting lo = lu-i +U0 in 
this equation, we find: 



(B.20) 



where S" is a term directly resulting from the frequency 
shift: 



S = 



r^-TB^ k^^Vr - {po{rVryy/po 
ikr [Por^-'Y 



kl + k-^ r-^^iBl 



(B.21) 
(B.22) 



The second line has been written for k^ >> 1/r. Because 
of S, the general compression depends on the background 
profiles. However, it can be noted that in the above 
limit {krT » 1), the frequency shift is very small and 
S l/{krr) << 1. Eventually, when neglecting the vari- 
ation of the background quantities on wavelength scale, 
but keeping the curvature effects, it is found that: 

V.v = 2— (B.23) 
r 

The Alfvenic perturbations do have a compression. The 
geometric interpretation is presented in the main body of 
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the paper and in fig. |21 The factor 2 comes from the fact 
that we require an initial equiUbrium state before perturb- 
ing with waves. The curvature of field lines exerts a force 
that must be balanced by a magnetic or thermal pressure 
gradient. This additional force is equal and opposite to the 
curvature force in the equilibrium state but acts together 
with the curvature in the wave equation to strengthen the 
compression. 

Remembering that the viscous dissipation does not de- 
pend exactly on the compression, we need to determine the 
quantity: 

D = V .V - 3d\\v\\ (B.24) 
From eq. IB. 81 we get the parallel contribution: 



im m? c? 



= —^0 = T^^V.^ (B.25) 
so that 

Z,= (l-3|)v.. = 2(l-3|)^I (B.26) 

From ea. lB.8l and lB.9l we see that the general Aflven waves 
keep their usual polarization: when fc^ = 0, if ve and Vz 
remain finite, then Vr must vanish, so that we always have: 

k^.vi_ ^ (B.27) 

As a consequence, when the Alfven waves are polarized 
perpendicular to the curvature direction, the curvature of 
the field lines has no effect on the compression and there 
is no associated viscous dissipation. However, the Galactic 
molecular clouds excite waves polarized almost isotropi- 
cally in every direction, so that the typical perturbed ra- 
dial velocity can be estimated of the order of the cloud 
speed Vc- Finally, defining a local radius of curvature Rc 
for the field lines closest to the cloud, we get in order of 
magnitude: 



« 2 ( 1 - 3 J- ) (B.28) 



At low-/3, D is dominated by the perpendicular divergence 
of the fiow whereas at high-/?, it is dominated by the par- 
allel contribution that may become very large for weak 
fields. When c^/v\ « 1/3, D vanishes and so does the 
dissipation. 



